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Abstract 

The thermal remagnetization (TR), i.e. the reentrance of magnetization upon heating in a 
steady-field demagnetized sample, is a common feature to the four types of polycrystalline per- 
manent magnets, mainly utilized for practical purposes, i.e. barium ferrites, SmCo§, Sm^Con 
and NdFeB magnets. The effect is small for pinning controlled and large for nucleation con- 
trolled magnets. The effect is strongly dependent on the demagnetization factor and may reach 
nearly 100% in SmCo^ samples measured in a closed circuit. The TR is very sensitive to a small 
superimposed steady field. The maximum effect and the position of the peak is dependent on 
the initial temperature. The direction of the TR is correlated with the temperature coefficient 
of the coercivity, resulting in a inverse TR in barium ferrite. The susceptibility of the thermally 
remagnetized samples is increased. Repeated cycles of steady-field demagnetization followed by 
heating result in the same TR. The phenomenology of TR and ITR is explained by means of a 
model taking into account both the internal field fluctuations due to grain interactions and the 
decay of single domain grains into multi-domain state. By taking the measured temperature 
dependencies of the coercivity and the saturation magnetization the theory is able to reproduce 
the experiments very well, allowing to determine the width of the field fluctuations, the width 
of the switching field distribution and an internal demagnetization factor as characteristics of 
the materials by fitting. 
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1 Introduction 

The effect of Thermal Remagnetization (TR) 
was first mentioned in a paper of Lifshits, 
Lileev and Menushenkov [I] devoted to stabil- 
ity investigations of the coercivity of SmCo§- 
magnets at elevated temperatures in 1974. 
They wrote: "We also found a remagnetiza- 
tion upon heating in absence of an external field 
at samples, which were partially or completely 
demagnetized by a steady counter field." 
One year later the same group [3] published 
a first systematic study devoted solely to this 
effect. As we will see in the following, TR is 
not restricted to SmCo^ magnets, instead it 
is quite common to all sintered hard magnets. 
But, before giving an outline of the history, we 
have to explain the principle of the TR exper- 
iment. In Fig. ^ the experimental procedure 
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Figure 1: The scheme of the standard TR- 
Experiment. After saturation the sample is de- 
magnetized isothermally by a steady field at 
the initial temperature To (upper sketch) and 
afterwards either heated (TR-experiment - left 
shaded picture) or cooled (ITR-experiment) 
with the external magnetic field switched off, 
whereby the magnetization is measured. 



is sketched. Firstly the sample has to be sat- 
urated by a very large magnetic field. Next 
a steady field Hr (remanence coercivity) op- 
posite to the saturation direction is applied. 
Switching off the external field drives the sam- 
ple in the dc-demagnetized state. Now the tem- 
perature is changed and the related change of 
magnetization is measured. If a reentrance of 
magnetism with increasing temperature is ob- 
served the effect is called thermal remagnetiza- 
tion (TR), see the right shaded picture. Other- 
wise, if a magnetization return is observed upon 
cooling, as shown in the lower shaded figure, we 
call it inverse thermal remagnetization (ITR). 
Already in their first investigations Lileev et 
al. PJ 13] demonstrated for SmCo^, that the 
TR is strongly dependent on the demagnetiza- 
tion factor, especially it amounts up to nearly 
100 % in a closed circuit. Furthermore it was 
shown, that at least for SmCo^, the effect is not 
necessary restricted to the demagnetized state, 
but also occurs in a partially dc-demagnetized 
state (± 23 % of the saturation magnetization). 
On the other hand the effect is absent in ac- 
demagnetized samples. Reconsidering the TR 
in SmCos nearly ten years later, Livingston 
and Martin f|] suggested as explanation of the 
TR the conversion of the single domain state 
(SDS) of some grains into a multidomain state 
(MDS). This assumption was supported by the 
observation of multidomain states in a few big 
grains of SmCo^ at higher temperature and 
also by the discovery of the TR in sintered 
AdFe-B-magnets [5|, where the effect is larger 
in samples with lower number of multidomain 
grains in the dc-demagnetized state. An alter- 
nate explanation was given in Ref. jSJ. Based 
on investigations of the TR in SmCo^ and bar- 
ium ferrite in 1985 the authors proposed the 
magnetic interaction between adjacent grains 
together with a broad distribution of the tem- 
perature dependent switching fields to be re- 
sponsible for the TR. This model, despite of 
its oversimplifications, was able to explain the 
increase of magnetization qualitatively and it 
showed the key role of the temperature coeffi- 
cient of the coercivity resulting in the discover- 
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ing of the ITR for sintered barium ferrites. Fur- 
thermore it was demonstrated that the TR is 
connected with a increasing susceptibility and 
that small external fields may suppress the ef- 
fect at all. Additional support the interaction 
model got from a paper of Lileev and Steiner 
0, proving that the TR effect depends strongly 
on the volume packing density, a fact which was 
shown to be responsible for the difference in 
ITR in sintered and chemical precipitated bar- 
ium ferrite samples 0. To decide, which of the 
both models accounts well for the observed ef- 
fects Lileev et al. [HI EH counted carefully grains 
in single-domain state and multi-domain state 
at different temperatures. Their results proved 
that there is a delicate interplay of both inter- 
action and formation of multi-domain grains, 
at least in the SmCo^ samples. From the 
very beginning it was assumed that TR is con- 
nected with the nucleation controlled magnets. 
Otherwise, pinning controlled magnets should 
not show the effect, and indeed, the SmCo^- 
, NdFeB-, and ferrite- magnets show remark- 
able TR and ITR resp. H 03 El EH HH EES , 
whereas the effect in Sm^Coyj is small [PH] . 
To what extend this statement is true together 
with the role of interaction in relation to the 
decay of SDS grains into MDS grains, moti- 
vated much of the theoretical and experimen- 
tal work. Since the effect is most conspicu- 
ous in SmCo*,, much effort was spent to this 
material. Otherwise it became clear, that it 
is a rather common effect to all modern poly- 
crystalline permanent magnets. Nevertheless, 
as will be shown in the next section, the ex- 
periments demonstrate clearly the individual- 
ity of every compound class. Regarding the nu- 
cleation controlled magnets we find that both 
the saturation magnetization and the coerciv- 
ity decrease monotonously for SmCo^- and 
NdFeB-magnets, but, while in NdFeB mag- 
nets the coercivity and the saturation magneti- 
zation vanish at nearly the same temperature, 
in SmCos the coercivity reaches zero at con- 
siderable lower temperature than the satura- 
tion magnetization. In contrast barium ferrite 
shows an increasing coercivity up to 550 K, i.e. 
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Table 1: Some of the main characteristics of the 
four groups of sintered permanent magnets. 

the temperature coefficient of the coercivity has 
the opposite sign. For the 5m2Coi7-magnets 
Ms(T) and Hq(T) behave as in S171C05, but 
there we have a complete different coercivity 
mechanism, i.e. volume pinning of domain 
walls. In NdFeB magnets the magnetic ac- 
tive phase is separated by one or more other 
phases reducing the interaction between adja- 
cent grains. In Tab. \I\we show schematically 
how the four groups of magnetic materials are 
qualitatively distinguished by some features be- 
ing relevant for TR, i.e. the number of phases, 
the coercivity mechanism, and the temperature 
coefficient of the coercivity. 

In the third chapter, we will give a brief re- 
view over the models used so far to understand 
the different aspects of TR and present a the- 
ory which accounts for both the TR in SmCo^ 
and the ITR in barium ferrite in single-phase 
magnets. The last chapter will be summarize 
what can be learned from TR measurements 
and some remarks to open questions will be 
appropriate. 

2 Experiments 

2.1 The standard TR experiment in 
the SmCo 5 , Sm 2 Coi 7 and NdFeB 

Sintered SmCos-magnets for which TR 
experiments are reported are single-phased 
and exhibit a poly-crystalline structure with 
the typical grain sizes much bigger than the 
critical radius for single domain behavior. 
They are mainly well textured. The phase 
diagram |141 I15j of the Sm-Co system shows 
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that below about 750 °C the magnetic active 
phase SmCos is meta-stable with respect to 
a decomposition into Sm,2Co-j and Sm^Con 
with reduced uniaxial anisotropy. While at 
moderate temperatures this decomposition 
process does not affect the microstructure and 
the magnetic quantities, chemical instability 
becomes a problem in high temperature mea- 
surements, since a long lasting TR experiment 
may act as an unintentional heat treatment 
of a sample. The first measurements of 
the TR in SmCo^ are collected in Fig. |2] 
The points shown as squares are taken from 
Ref. [3] and the triangles from Ref. [I] and 
the open circles are our own measurements 
at a VACOMAX 200 specimen [IE]. The 
biggest TR (filled squares) was measured in a 
closed magnetic circuit what corresponds to 
zero demagnetization factor. The remaining 
curves represent open circuit measurements 
at cylindrical samples of different aspect 
ratio resulting in different demagnetization 
factors. In Ref. [3] the remanence was mea- 
sured at room temperature after heating to 
a certain temperature and cooling down to 
room temperature, whereas in Ref. 0] the 
magnetization was registered simultaneously 
while increasing the temperature stepwise. 
Every temperature step was followed by an 
annealing of 30 min. The remaining curves 
were measured while continuously increasing 
the temperature. In order to get the slightly 
differing experiments into one diagram we used 
the temperature dependence of Ms measured 
at our VACOMAX 200 sample together with 
the room temperature remanence data given in 
Refs. |1] to scale the results. Although one 
may object that the results of Fig. Ostein from 
SmCos samples produced in different ways, it 
is evident that a variation of the demagneti- 
zation factor influences the TR much stronger 
than a variation of the production conditions, 
at least if the result is a good SmCo^ hard 
magnet. Fig. |2] shows both the TR for two 
samples with different demagnetization factors 
cut from the same material (VACOMAX 170) 
[Tfij . and the temperature dependence of the 
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Figure 2: The TR for various SmCo^ sam- 
ples taken from Ref. [3] (squares) and Ref. [1] 
(triangles) and a VACOMAX 200 sample (cir- 
cle) [Ej- The data were transformed slightly 
from the original data to fit into the same plot. 
Also the demagnetization factors of the cylin- 
drical samples were calculated from the geom- 
etry given by the authors. For details we refer 
to the original papers [31 Q]- 

saturation magnetization and the coercivity of 
that material. It can be seen from the curve 
for the N=0.08 sample that the TR survives in 
the region above the temperature Th c , where 
the coercivity vanishes, which is about 700 K. 
As we will see in the theory, this unexpected 
behavior can be ascribed to the fact, that the 
Curie temperature is more than 200 K higher 
than Thc- Lileev et. al. measured the TR 
for SmCo§ samples with coercivities between 
(0.35 T/7/o and 3 T/fiQ. They found that the 
inflection point, visible also in Fig. |3] vanishes 
for room temperature coercivities J3| below 
1.6 T/hq. For higher coercivities they found 
a decrease of the initial slope and an increase 
of the maximum of the TR with increasing 
coercivity. Measurements at samples, where 
the low-temperature coercivity was reduced 
by a special heat treatment ^H] but otherwise 
prepared in the same way as the samples 
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Figure 3: The temperature dependence of the 
coercivity He, saturation magnetization Ms 
and the TR for two samples of SmCo^ with de- 
magnetization factors N=0.08 and N=0.53 [T5| . 




100 200 300 400 500 600 700 800 900 



T/K 

Figure 4: The temperature dependence of the 
coercivity He, saturation magnetization Ms 
and the TR for a SmiCoyj sample. 



utilized in Fig. |31 did not show significant 
changes of the TR curves. 

In contrast to SmCo^, where the hystere- 
sis mechanism is nucleation controlled, the 
coercivity of S1712C017 magnets is caused by 
volume pinning of domain walls. For that be- 
have a single phase solution of approximately 
1:7 composition, i.e. slightly Co deficient 
with respect to SmzCon, is homogenized at 
approximately 1200 °C and afterwards cooled 
down to 800 °C, where it decays into the 
two adjacent stable phases by precipitation 
of SmCoz and SmCu^. The contemporary 
commercial magnets of the 2:17 type are 
composed typically Sm(Co,Fe,Cu,M)7, where 
the iron is used to increase the saturation 
magnetization and the other metals mainly for 
metallurgical reasons to obtain the precipita- 
tion microstructure. The latter is common to 
all £m2Coi7-magnets. The entirely different 
mechanism of magnetic hardening was the 
reason that in early years these compounds 
where believed not to show TR. Nevertheless 
a small TR is present also in SV712C017, as it is 



visible in Fig. 0] With regard to the smallness 
of the effect, which amounts to not more than 
3%, it is not yet clear, whether the effect is 
characteristic for all Sm^Coyj magnets or 
whether it is caused by imperfections of the 
microstructure, especially by the presence of 
a small content of of SmCo^ phase, as was 
argued in Ref. ^Hj m connection with the 
interpretation of virgin curves. 

The NdFeB-type magnets consist of the 
magnetic active phase Nd2- x Dy x Fe\4 : B sur- 
rounded by nonmagnetic phases. This guar- 
antees for magnetic decoupling of neighboring 
grains. Within the main phase the nucleation 
mechanism was shown to control magnetization 
reversal |19j . Apart from the much lower Curie 
temperature and the spin reorientation below 
137 K, the temperature dependence of coer- 
civity and saturation magnetization is qualita- 
tively the same as in SmCo^ within the inter- 
esting temperature region, see Fig. |3 A re- 
markable difference is that coercivity and mag- 
netization vanish nearly at the same temper- 
ature, i.e. at T}ic~ Te- The observed TR, 
as depicted in Fig. 03 is considerably smaller 
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Figure 5: The temperature dependence of 
the coercivity He, saturation magnetization 
Ms and the ITR of two NdFeB samples with 
N=0.04 and N=0.46 resp. QS]. 

than in SmCo^, but, with 12 % for the sam- 
ple with very small demagnetization factor it 
is doubtless that it is an intrinsic character- 
istic of the magnet and can not be ascribed 
to imperfections in the microstructure. Liv- 
ingston reported TR effects of 3% for rod sam- 
ples (N»0.17) with H c about 1 T//i and 24 % 
for a high coercive sample (2.5 T//xo) 0- Fur- 
thermore, the low coercive samples were shown 
to contain a lot of multidomain grains in the 
dc-demagnetized state, in contrast to the high 
coercive specimen, which contained only a few. 
The same correlation between coercivity and 
TR was also reported in [TT] . 

2.2 The ITR in Barium Ferrite 

In comparison with the above mentioned in- 
termetallic magnets the ceramic barium ferrite 
magnets exhibit a entirely different tempera- 
ture dependence of the coercivity, as shown 
in Fig. El That is why the normal TR can- 
not be measured in this substances. Nev- 
ertheless an analogous experiment, where an 
initially saturated sample is dc-demagnetized 




T/K 

Figure 6: The temperature dependence of the 
coercivity He, saturation magnetization Ms 
and the ITR of a barium ferrite sample with 
N=0.1 

at elevated temperature and afterwards cooled 
down, shows also a remarkable remagnetization 
jSJ. This effect was called inverse thermal re- 
magnetization (ITR). As seen from Fig. Qwc 
find again a strong dependence on the demag- 
netization factor. Since the coercivity upon 
cooling does not vanish one cannot observe 
a maximum. Instead one finds a plateau at 
lower temperature. Barium ferrite magnets are 
single-phased and the coercivity mechanism of 
barium ferrite is nucleation controlled. In that 
respect, apart from the different temperature 
dependence of the coercivity, the ITR of about 
20% is similar to the TR in SmCos-magnets. 

2.3 The field dependence 

The dramatic dependence of the TR on the 
demagnetization factor shows that a high sen- 
sitivity with regard to the internal magnetic 
fields is to expect. The influence of a small su- 
perimposed steady field was studied in [HJ 12*01 
12 lj . To that behave the standard experiment 
was modified in the following way: After sat- 
urating the sample in a field of (about 10 T 
for SmCo^- and 5m2Coi7-samples and about 
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Figure 7: The ITR for four samples with 
N=0.1, N=0.3, N=0.7, and N=0.9 respectively. 
The solid lines are theoretical fits. 



2 T for the barium ferrite samples) it was de- 
magnetized by help of an opposite steady field 
Heirt- -^ this field equals the remanence coer- 
civity Hr the magnetization will go to zero af- 
ter switching off the field. Otherwise, if H^ x \ 
deviates a little bit from Hr, the demagne- 
tized state will be achieved on the recoil curve 
for a small negative residual field H ext < 0, 
if H^ xt is greater than Hr, and for a small 
positive field otherwise, as the inset in Fig. |H1 
shows schematically. The field H ex t was then 
kept constant while the sample was heated (or 
cooled in the case of barium ferrite) and the re- 
sulting remanence enhancement AM(T) (TR- 
curve) was recorded. A set of such measured 
TR-curves (points) is given in Fig. [S|for SmCos 
[221 ■ Tne initial temperature T was 250 K. The 
solid lines are due to the theory, which will be 
explained below. The cases with H ex t = 0, 
i.e. Hi = —Hr, correspond to a "normal" 
TR-curve. A systematic shift of the maximum 
temperature T max , the maximum remanence 
enhancement AM max and also of the initial 
slope AMtr(T)/ AT 22 with increasing nega- 
tive external field is observed. The same qual- 
itative behavior can be measured in NdFeB- 
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Figure 8: Measured TR-curves AM(T) (cir- 
cles) for SmCos |22| for different external 
steady fields applied during heating (small 
numbers at the curves). The solid lines are 
model calculations with parameters determined 
by the H ex t = 0-curve (cf. section 3). The in- 
set shows schematically how the initial state for 
the TR measurement is prepared. 



magnets (Fig. EJ) and SraiCoyj magnets (Fig. 
I10|) . Besides the fact, that the absolute TR- 
change with field is smaller in NdFeB and 
even in Sm^Coiy-samples, it looks very simi- 
lar if the remanence-change with temperature 
is related to the remanence enhancement of the 
standard experiment ^) . Fig. shows the 
analogous experiment for the ITR measured at 
a barium ferrite sample [22] . As it is obvious 
from the Figs. 01101 (for NdFeB cf. [22]) a 
small positive field increases the TR extremely. 
Otherwise the TR may be suppressed at a cer- 
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the curves) [50] . 



(small numbers at the curves). 



tain negative field H' . The order of magnitude 
of this "suppression field" H' is to some ex- 
tent a characteristic of the different types of 
permanent magnets [0J [50] . Fig. ^] extracts 
the maximum TR values in dependence of H ex t 
from Figs. [§] [0] ^] and fTTl Comparing these 
results with other reported measurements one 
finds H' of the order of 200-500 mT for SmCo 5 
H0J HI] and 30-80 mT for NdFeB -magnets 
[2UH2]. For Srri2(Co, . . .)i7-magnets this value 
seems to be more dependent on the chemical 
composition and technological treatment of the 
sample. Please note the mT-scale in the lower 
picture of Fig. ^] showing AM max via H ext 
data for Sm^Coxj. Nevertheless H' is in the 
same region as in SmCo^. This may support 



the assumption that the TR in Sm^Coyj is 
caused by a small amount of SmCo^ grains as 
was proposed in [T^]. The suppression field for 
the ITR in the barium ferrite sample of Fig. ^2 
is about 20 mT as may be seen from the filled 
squares in Fig. 1121 



2.4 Influence of the initial tempera- 
ture T 

Ocularly the TR as well as the ITR are strongly 
connected with the temperature dependence of 
the coercivity. The difference of the coercivity 
Hc(T max ) — He (To) is essentially for the TR. 
Thus the remanence increase should be depen- 
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Figure 11: Measured ITR-curves AM(T) (cir- 
cles) for barium ferrite for different exter- 
nal steady fields applied during cooling (small 
numbers at the curves) 



dent on the initial temperature Tq as well. In 
Fig. El we show the results for oblate discs 
cut from a well aligned SmCc^-magnet (VA- 
COMAX200). The upper part shows the initial 
irreversible TR at very low temperature, i.e. in 
the high coercive region, in detail, whereas the 
lower part shows the complete measurement to 
the maximum. Cooling down from the maxi- 
mum the magnetization increases reversibly - 
exactly proportional to Ms{T). More gener- 
ally, if the TR experiment is interrupted at 
some temperature T\ the magnetization change 
while cooling follows exactly mi x Ms(T), with 
mi being the ratio AM(Ti)/Ms(Ti), see e.g. 
Fig. 1 in 20 . On the one hand the amount of 
the TR is considerably influenced by a change 



0.5 

0.4 
H 0.3 h 
J 0-2 
=*- 0.1 


-0.1 



i i i i i i i i i i i i i i i i i i i i i i i i i ; i i i 




I i i i i I i i i i I i i i i I i i i i I 



-250 -200 -150 -100 -50 50 
u.H f /mT 

~0 ext 



s 




-250 -200 -150 -100 -50 50 



"0 e 



/mT 



Figure 12: The maximum TR-curves in depen- 
dence of the superimposed magnetic field for 
the samples of Figs. |H1 El El EH Open trian- 
gles: SmCo5, filled triangles: SmiCoyj, dia- 
monds: NdFeB, squares: Barium Ferrite. 

of the initial temperature, on the other hand 
only a small shift of the maximum position 
to lower temperatures is observed, if the ini- 
tial temperature decreases. A likewise behav- 
ior is observed in NdFeB as depicted in Fig. 
1141 Regarding the ITR in barium ferrite, where 
the initial temperature has to be chosen higher 
than the final temperature, we found again a 
strong dependence of the remanence increase 
on the initial temperature, whereas the region 
where the curves run into the plateau is slightly 
influenced only, as may be seen from Fig. 1151 

2.5 Susceptibility and Repeating ex- 
periments 

Since heat treatment is a mayor step in 
permanent-magnet production the possibility 
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Figure 13: TR curves for different initial tem- 
peratures To. The upper picture shows the ini- 
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is to envisage that the sample after the first 
TR experiment is irreversible changed. This is 
of course of minor importance for the ceramic 
magnets but may be important in SmCo^. 
Since the maximum temperature in SmCo§ is 
about 600 K, it is low enough not to destroy the 
samples if the TR experiment does not last to 
long. In [3] the stability of the magnetization of 
SmCos against thermal cycling between 20 °C 
and 75 °C was investigated. It was shown that 
it needs a few cycles to stabilize the magnetiza- 
tion change. For a sample which was partially 
dc-demagnetized (10% lowered remanence), it 
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Figure 14: TR curves for two different initial 
temperatures for a NdFeB sample (NdFeB-HR 
510) with N=0.46 



was shown that 16 h annealing at 75 °C after 
ten cycles produces a sharp positive AM jump 
while additional hold of 38 h after cycle 13 pro- 
duced no further significant change, besides a 
very small increase continuing up to 53 cycles. 
In [22] "repeating" experiments were reported, 
which differ a little bit from the thermal cy- 
cling procedure. Contrary to [1], where after 
every TR cycle the sample was saturated with 
a high field, the sample was dc-demagnetized 
from the remanent state which it reached after 
a cycle To — > T max — > To, i.e. without a new 
saturation. The idea behind this experiment 
was to prove the assumption, that the "hard 
grains" , i.e. the fraction of grains which resist 
dc-demagnetization, are not concerned by the 
TR experiment. Fig. ED shows a TR exper- 
iment together with its two repetitions. The 
difference between the curves are more likely 
due to the uncertainty in achievement of the 
dc-demagnetized state. Thus it is merely the 
"weak" grain fraction, which is switched both 
during dc-demagnetization and TR. After the 
first heating the sample is demagnetized by a 
much lower field, as may be seen from the in- 
set of Fig I16| showing that the magnetic hard- 
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Figure 15: ITR curves for three initial temper- 
atures for a isotropic sample. 

ness of the week grains was considerably low- 
ered by heating to T max . As a consequence an 
increase of the susceptibility is observed. In 
the experiment shown in Fig. El the TR curve 
was started as normal in the dc-demagnetized 
state, but a very small field of 10 mT was su- 
perimposed with alternating sign at every reg- 
istration point of the magnetization. This re- 
sults in a "normal" and two slightly shifted 
TR curves. The latter were taken as input for 
the calculation of the susceptibility x, which 
was derived from the un-sheared susceptibility 
X ' = AM/AH ext according to X = x'/(l-iVx') 
with TV being the demagnetization factor of the 
sample. The plot in Fig. 1181 shows a continu- 
ous uprise with increasing temperature, which 
becomes very large beyond T max . This sharp 
increase is due to the partial conversion of the 
weak grain fraction into the MDS and due to 
the lowering of their switching fields. This re- 
duced coercivity of the weak grain fraction is 
frozen if the sample is cooled back to To. Oth- 
erwise, the hard-grain fraction regains nearly 
reversibly its initial state, what means that the 
interaction field distribution due to the hard 
grains is-reestablished. In order to demagne- 
tize the sample another time one has to over- 
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Figure 16: TR curve (1) and two repetitions 
(2,3) for the same sample of SmCo^. The in- 
set shows the dc-demagnetization curve for the 
first and second run. 

come the smaller coercivity of the weak-grain 
fraction only. Afterwards the switched weak 
grains are located in nearly the same environ- 
ment as before. As a consequence the TR is 
nearly the same in the repeating experiments. 

2.6 Influence of texture. 

Due to the technical importance, most of the 
reported TR experiments were done at well- 
textured samples. In Refs. [221 ESI TR in- 
vestigation at isotropic magnets were reported. 
Comparing isotropic and well aligned samples 
of SmCo^ one finds the same qualitative be- 
haviour, but the TR of the isotropic specimen 
is reduced roughly four times with respect to 
the well aligned sample. Furthermore the dif- 
ference of the TR curve and the first repeti- 
tion differs a little bit more than in the well 
aligned material, as may be seen from Fig. ED 
in comparison with Figs. 12131 and ^3 Qual- 
itatively the same behavior was reported for 
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Figure 17: Measured TR curves for the same 
sample as in Fig. [HI for fiQH ext = ±10mT (+). 
The circles represent the values measured for 
fioH ext = 0. Furthermore the Hq{T) is indi- 
cated. 

isotropic sintered barium ferrite Otherwise, 
an isotropic barium ferrite sample made from 
chemical precipitated 50 nm-powder, does not 
show any TR [HJ To what extent this 

is connected with texture is not completely 
clear. The nano-sized grains are nearly Stoner- 
Wohlfarth particles, with an entire different co- 
ercivity mechanism (2Z1- Furthermore the den- 
sity is much lower than in the sintered barium 
ferrite samples, what reduces the inter-grain in- 
teraction. That the latter plays an essential 
role can also be concluded from the experi- 
ments of Lileev and Steiner [7j , who measured 
the initial TR for SmCo^ in the temperature 
region from 4.2 K to 270 K for both well aligned 
sintered magnets and polymer-bounded sam- 
ples with a SmCoz, content of 20% and 80%. 
They found a considerable reduction of the TR 
with decreasing SmCo§ packing density. 

3 Theory 

3.1 Early models 

During the TR experiments the macroscopic 
magnetic energy is increased. This is only pos- 
sible if the internal energy is considerably de- 




T/K 



Figure 18: Measured and calculated suscep- 
tibility (solid line) fj,oH ex t = ±10mT for the 
mCos-magnet from Fig. El Furthermore the 
theoretically calculated percentage of "weak" 
grains is shown. 



creased at the same time. During the initial 
dc-demagnetization process a state is prepared 
where one half of the grains is upwards (i.e. in 
the initial saturation direction) and the other 
half is downwards magnetized. Under the as- 
sumption that the intrinsic switching fields of 
the grains are not strong correlated with its po- 
sition, we will find a extreme inhomogeneous 
magnetization at a microscopic scale in the dc- 
demagnetized state giving rise to considerable 
stray fields. While there is agreement that it is 
this part of the energy which governs TR, two 
different explanations for the mechanism be- 
hind TR were proposed. The first, introduced 
by Livingston [1], assumes that the downwards 
magnetized grains, which are in a single domain 
state before heating, decay into a multido- 
main state upon heating. This model got sup- 
port from the direct observation of such grains 
in Kerr photographs taken from the surface 
of dc-demagnetized S'mCos-magnets [4 and 
iV(iFe.B-magnets 0. Although these intrigu- 
ing observations make clear that to some extent 
nucleation of domain walls are involved in the 
TR, there are also some arguments against the 
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Figure 19: TR curve (1) and one repeti- 
tion (2) for an isotropic sintered sample of 
SmCo 5 (Quality as VACOMAX 170, sintered 
and pressed without aligning field |IB])- The 
inset shows the dc-demagnetization curve for 
the first and second run. 



dominating role of the SDS-MDS mechanism. 
Without any calculations it is obvious that 
the maximum TR effect may not be greater 
than 50%, since after a complete decay of the 
weaker half only the harder grains contribute 
to the magnetization. This contradicts the ob- 
servation of a much greater effect measured 
for SmCos in a closed circuit [3]. Further- 
more, the number of the multidomain grains 
seen in the shown Kerr photographs seems to 
be to low to explain the effect, despite the 
fact that it is not completely clear to whether 
extent the surface observation is representa- 
tive for the bulk behavior. The other model 
[H] assumed that the grain fraction with lower 
switching fields is first switched downward by 
the dc-demagnetization procedure. After tak- 
ing away the external field they remain in a 
metastable state since their switching field is 
higher than the interaction field caused by the 
surrounding upwards oriented, i.e. harder mag- 
netic, grains. Upon heating the switching fields 
are reduced and the weaker grains are switched 



back in its initial direction. The decrease of 
the TR with decreasing filling factor [Jj hints 
strongly at the role of inter-grain interaction. 
Further support the interaction model got from 
the ITR in barium ferrite, since contrary to 
SmCo^, S1TI2C017, and NdFeB, the low tem- 
perature minimum of the coercivity is still high 
enough to prevent a decay of grains into a 
multidomain state, thus the SDS-MDS mech- 
anism can not explain the ITR. To clarify the 
role of multidomain grains Zaytzev et. al. [2j 
used Kerr micrography to count the number 
of multidomain grains in SmCo^ magnets in 
the dc-demagnetized state at room tempera- 
ture and at 120 °C. They demonstrated clearly 
that in the considered samples three processes 
contribute to the TR: Downward magnetized 
grains are either switched back or decay into 
the MDS under the influence of the interaction 
fields due to surrounding grains and in MDS 
grains the upward domains grow at cost of the 
downward domains. Besides the fact that the 
investigated samples exhibit a lowered coerciv- 
ity and an unusual broad switching field distri- 
bution (SFD) in comparison with usual com- 
mercial magnets and that surface effects may 
result in a slight overestimation of the role of 
the MDS fraction from these experiments it be- 
comes clear that a theoretical model has to de- 
scribe both the interaction processes and the 
behavior of the MDS grains. It is interesting 
that they also observed hard grains changing 
from SDS to MDS. Whereas the theory pre- 
sented below uses a phenomenological ansatz, 
Lileev et al. |T7] presented a model based on 
simulations of interacting dipoles. This model, 
also taking into account the possibility of MDS 
states, was originally developed for the simula- 
tion of minor loops of permanent magnets [2E] • 
By taking into account the temperature depen- 
dence of the model parameters it was also pos- 
sible to simulate the TR. The shape of the only 
given curve deviates far from experiments. The 
authors see the reason for that in the mean- 
field treatment of the interaction contribution 
within their model. Whereas Monte-Carlo sim- 
ulations allow to play with parameters it is al- 
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most impossible to determine parameters by 
fitting to the experiment. For that purpose a 
phenomenological description as given in the 
following seems more appropriate. 

3.2 The inclusion model 

The outlined theory was developed in Refs. 
[El EE EH1 EH] • The poly-crystalline permanent 
magnet is an ensemble of high-coercive mag- 
netic uniaxial grains. To fix our model we have 
to define the properties of both a single grain 
and the ensemble. For the single grain we as- 
sume that the temperature dependent satura- 
tion magnetization M${T) should be known. 
Furthermore every grain should be character- 
ized by its temperature dependent switching 
field H S (T), by its "internal" demagnetization 
factor re, by the direction of the easy axis due 
to the high uniaxial anisotropy c, and by its 
volume V. Here H S (T) is the absolute value 
of the field needed for reversing the magnetiza- 
tion of the grain in a closed circuit. Regarding 
the ensemble we assume that all grains exhibit 
the same Ms(T), that the easy axis c of the 
grains are completely aligned, i.e. we restrict 
the model to ideally textured, single-phased 
magnets. These conditions are best fulfilled in 
SmCos and barium ferrite magnets, whereas 
NdFeB magnets contain more than one phase 
and Sm2Con magnets show a precipitation 
structure resulting in a different switching be- 
havior of a grain. Furthermore we allow that 
the grains differ in their switching fields H s , but 
the distribution of these switching fields should 
be known. In non-interacting ensembles the 
switching field distribution may be determined 
from remanence measurements [213] , but if the 
grains interact strongly this method fails. In 
Ref. [23] both the thermal remagnetization and 
the change of the H s spectrum with increas- 
ing temperature of a SmCo^ sample was mea- 
sured, whereby the latter was approximately 
determined from dMj dH. This is justified if re- 
versible magnetization processes are negligible 
as it is the case in SrreCos-magnets, For sim- 
plicity we assume a Gaussian normalized with 



respect to the region < H s < oo 10 , 
2 



v^fl + erf (f 



(Hs — Hs 



1) 



with the "mean switching field" H s and the dis- 
tribution width a s . The temperature depen- 
dence (not the absolute value!) of the switch- 
ing fields of all grains are equal. For the mo- 
ment we have to regard H s and a s as model 
parameters. Fortunately we will see later on 
that the theory delivers a possibility to relate 
H s to the coercivity Hq of the sample. An cen- 
tral point of TR is the treatment of the grain- 
grain interactions resulting in a collective be- 
havior of adjacent grains. On a macroscopic 
length scale the internal magnetic fields and the 
magnetization of the sample are homogeneous. 
But, if the length scale is reduced to the order 
of some grain diameters the magnetization be- 
comes coarser and the inhomogeneities gain in- 
fluence. Averaging over volume elements con- 
taining a few grains only, yields values, which 
deviate from the average. Thus, we consider 
every grain as an inclusion embedded in a local 
environment, which may differ stochastically in 
its magnetic field and magnetization from the 
related averaged values. We show this schemat- 
ically in Fig. 1201 For simplicity we approxi- 
mate the magnetization in the vicinity of the 
inclusion by the mean magnetization (M), ne- 
glecting the mentioned magnetization fluctua- 
tions. Otherwise we allow for field fluctuations 
AH = H — (H) around the mean internal field 
in the environment of a grain. The mean inter- 
nal field is related to the external applied field 
H ext according to 



f i (H) = fM) H ext -N(M), 



(2) 



with N being the demagnetization factor of the 
sample. The field fluctuations are assumed to 
be Gauss distributed [TU1 13*1] . 



f(H) 



1 - ¥f 



with 



(3) 



AH = H - (H) . 



13 



Inclusion (grain) 
Hi^H-nC-Ms -<M)) 
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Figure 20: The grain is considered as an in- 
clusion with magnetization ±M$ (here we as- 
sumed —Ms, symbolized by the long white ar- 
row) embedded in an environment differing in 
its magnetic field H from the mean magnetfield 
(H) within the sample due to the fluctuation 
AH. Fluctuations of the magnetization in the 
environment are neglected, thus the magneti- 
zation in the environment is equal to the mean 
magnetization of the sample. This is symbol- 
ized by the two shorter white arrows of equal 
length. 

what is a consequence of the central limit the- 
orem of probability theory. Furthermore, the 
fluctuation width &f has to be an even func- 
tion of the magnetization itself. 

(M) 



m 



with 



m 



M s 



(4) 



This becomes obvious since in the saturated 
magnetized state the fluctuation width is zero 
(remember: ideal texture, all grains the same 
phase). Otherwise it is maximum for the dc- 
demagnetized state with random distribution 
of upward and downward magnetized grains. 
A Taylor expansion of oy((M)) to second or- 
der which fulfils these two boundary conditions 
yields 

,2\ 



£( m ) = (1 - m z ) . (5) 

In general £(m) depends on the correlation of 
neighboring grains. For instance the field of 



a simple cubic array of uncorrelated dipoles, 
randomly oriented in ±z-direction, yields 



£( m ) = Vl 



(6) 



Of course the local fields fluctuate around the 
mean internal field also in direction, thus, 
strictly speaking eq. © accounts for the fluc- 
tuations of the z-component only. Furthermore 
the deviation of the field direction from the 
easy axis gives rise to rotation processes, but 
this is surely a small effect in ideally aligned 
magnets with high anisotropy constants. In 
the early works on the TR theory a further 
point was stressed (HJEIlEh The field depen- 
dence of the switching field. This is indeed a 
crucial point, if one attempts to describe re- 
peating experiments or the virgin curve, but 
can be neglected in the first TR experiment, if 
Hji is large enough to make most of the down- 
switched particles hard, as it seems to be the 
case in SmCo§ and barium ferrite. Taking into 
account this process is in principle not difficult, 
since it simply scales the SFD of the switched 
particles. This would introduce at least one ad- 
ditional parameter. We abstained from doing 
so by restricting the theory to the first TR run 
only. 

3.3 Calculation of the Isothermal 
Demagnetization Curve 

The first step is the calculation of the demagne- 
tization curve from the saturated state down to 
the magnetic field H^ t in the third quadrant, 
see Fig. ^ To save the opportunity to de- 
scribe also the influence of a small steady field 
or a non-zero initial remanence on the TR we 
use a more general scheme of the TR experi- 
ment as shown in Fig. 1211 After the saturation 
all grains are up (f) magnetized. They posses 
their maximum switching fields and (M) = Ms 
holds. If the external field is lowered to a value 
^ext ( c -f- Fig- [^the magnetization of a grain 
is switched down (J.) if its internal field Hf) 
is lower than — H s and the internal field after 
switching H^) will be smaller than +H S . Thus 
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Figure 21: Generalized scheme of the a TR- 
Experiment: After saturation the sample is de- 
magnetized isothermally by a steady field to an 
arbitrary point of the demagnetization curve, 
indicated by H\ and M\. Afterwards the field 
is reduced to an arbitrary value Hi- The re- 
lated magnetization is M2. Finally the sam- 
ple is heated (TR-Experiment) or cooled (ITR- 
Experiment) with the external magnetic field 
kept constant. 



we have the two conditions 



H 



(i) 

i,1 



< 



and 



H^<+H s .(7) Pw{Hs 



The internal fields before and after switching 
are with respect to the above introduced ap- 
proximations 



H 



(i) 



,,, F-n(M s -(M)x) 
and 



H 



(i) 

hi 



H — n(—M s — (M)x) 



(8) 



(9) 



Here H is the local field in the environment 
of the grain. A variation of the external field 
H ex t results in a change of (M) on the one hand 
and on the other hand it changes the probabil- 
ity f(H)dH that the field in the environment 
of a grain is between H and H + dH due to 
the relation © and eq. ijljl. If the first condi- 
tion from eq. (JZJ) is fulfilled but not the second, 
the grain cannot jump into a stable state. If the 
grain is large enough it solves the conflict by an 
incomplete jump, what turns its single-domain 



state (SDS) into a multi-domain-state (MDS). 
We will call such grains "weak". Whereas the 
"hard" grains can exist in states with ±M$ 
only, the magnetization of a weak grain (M)i 
is an average over the upwards and downwards 
magnetized volume fractions within the grain. 
From phase theory of 180°-domains we get for 
the averaged magnetization of the i th grain in 
dependence on the local field in the environ- 
ment 



(M) + - 
n 



(10) 



Opposite to the hard grains the weak grains 
have no memory since the magnetization fol- 
lows the local field immediately as long as 
|(M)j| < Ms holds. Otherwise they are sat- 
urated up- or downwards. If one calculates 
the probability that a grain with a given H s 
is in a MDS the switching condition (J7| limits 
the integration over the local field distribution 
from above, whereas the back-switching con- 
dition H^) > +H S limits from below. The 
related probability is 



dH {(H) 



Hl 



with 

H H 
H L 



-H s + nM s - n(M) 
H s - nM s - n(M) . 



(11) 



(12) 
(13) 



The requirement, that the upper limit has to 
be greater then the lower one, yields the con- 
dition H s < nM s . Thus, grains with switching 
fields smaller than their own internal demagne- 
tizing field are weak. The magnetization of the 
volume fraction of the weak grains with a given 
H s is 

00 

(M) W (H S ) = M s J dHi(H) 
1 

j dHf(H) 



-M S 
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+ j dHi{H) (M)i 

Hr 



(14) 



The integral may be easy evaluated, where we 
have to take for the magnetic field H^ t and the 
related magnetization (M)i at the hysteresis 
curve. For hard grains Hl is higher than Hjj- 
If these grains fulfill the switching condition, 
the second condition on eq will be fulfilled 
automatically. The probability to find a grain 
downwards magnetized if a field H^) t is applied 
is 

+00 

Pi(H s ) = jdHf (1) (H) x 

—00 

xO(-H s - H$) &{H S - H$) , (15) 

with H5V and H?) being the fields within the 
grain in the upwards and downwards magne- 
tized state resp. and 0(x) is the Heaviside 
function, fm is the field distribution eq. (jlj 

with (H) = H^J t - N(M) 1 . After carrying out 
the integration in eq. (|15j) we get for the mag- 
netization (M)i = M s (1 - 2p x (H s )) of a frac- 
tion of hard grains with a given H a 



{M)l(Hg 

with 
x H = 



-M s eri(x H ) 



(16) 



(17) 



-H a + nM s - n(M} 1 - H^ t + N(M) 1 
it/ (mi) 

The total magnetization of the sample one gets 
from averaging with respect to H s 

nMs 

(M)x = J dH s g(H s )(M) wA (H s ) 



00 

+ j dH s g(H s ) (M)x(H s ) (18) 

nMs 

Here the integration has been split due to the 
different contributions of the weak and hard 
grains. For the weak grains we have to use eq. 



(|TH . Since the right hand side of eq. Q18|) de- 
pends on the mean magnetization (M)i itself, 
we have to solve this implicit equation numer- 
ically. 

3.4 Calculation of the Isothermal 
Recoil Curve 

Next, the external field is changed from H^ t 
to H^} t with H^) t > H^ t via the recoil curve. 
Some of the hard "down" -grains are switched 
back, whereas the weak grains shift their mag- 
netization reversible. The related mean mag- 
netization is (M)2- The contribution of the 

weak grains results from eq. (|14j) . if we in- 

(2) 

sert now H y e J t for H ext and (M)2 for (M) resp. 
followed by an integration from zero to nMs 
with respect to H s . The contribution of the 
hard grains is more difficult to calculate, due to 
the memory effect. Let us consider the proba- 
bility ioit(iT s ) that a hard grain with a given 

H s was switched by H^ t and switched back 

(2) 

by i^ext afterwards. To calculate this probabil- 
ity it is inevitable to make an assumption on 
the correlation between the fluctuation of the 
local field Hi which acts on a grain if H^) t is 
applied and the fluctuation of H2 according to 
the external field H^} t . Of course, if H^} t is 
only slightly different from H^ t it is unlikely 
that the neighborhood of a grain changes con- 
siderably, so that Hi and H2 should be strongly 
correlated. With increasing distance between 
H^ t and H^) t this correlation will vanish, due 
to the multitude of switching processes, hence 
we can average with respect to Hi and H2 inde- 
pendently, if the difference H^) t — H^ t is not to 
small. Thus the probability decouples ac- 
cordingly wii(H s ) = pi(H s )q 2 (H s ) withpi(iT s ) 
from eq. H15|) and 



Q2(H S ) 



DC 

/ 



dHf (2) (H) x 



x e(#g } - H a ) Q(H% + H s 



(2) 



1 - erf(y L ; 



(19) 
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with 
VL -- 



(20) 



H s - nM s - n(M) 2 - ffg) + N ( M ) 
&f(m) 



Due to the mentioned magnetization changes in 
the environment of a grain it may happen that 
grains which resisted H^) t are switched down 

(2) 

by H e J t ("up" -"down" -contribution). The re- 
lated probability factors also 

w U {H s ) = (l-pi(H s ))p 2 (H s ) (21) 

with p\(H s ) again from eq. (fT3|) and 

oo 

P2 (H S ) = J dHi {2) (H)x 



r(2)> 



with 
z H ■■ 



(22) 
(23) 



-H a + nM s - n(M) 2 - + N(M) 2 
&f(m) 

There are two further probabilities correspond- 
ing to the remaining two histories, which a 
grain may experience, i.e. the probability that 
it resisted both and H^ t ( "up" -"up" con- 
tribution) and the probability that a grain was 
switched by H^ t but resisted back-switching 
by H^ t ("down" -"down" contribution). We 
find 



and 



1- P1 (H S )\ [l-p 2 (H s )\ (24) 



P1 (H S ) (l - q 2 (H s 



(25) 



Thus hard grains with a given H s contribute to 
the magnetization 



M 2 {H S ) = M s (w n (H s )-w n (H s ) 



(26) 



Averaging with respect to H s yields the re- 
coil curve for the mean magnetization {M) 2 in 
dependence on both magnetic fields H^l and 



rr(2). 
n exf 



(M) 2 = J dH s g(H s )(M) Wt2 (H s ) (27) 
o 

oo 

+ J dH s g(H s )(M) 2 (H s ). 

nM s 

3.5 Calculation of the (Inverse) 
Thermal Remagnetization 

The above calculation of the demagnetization 
and the recoil curves was done with the implicit 
understanding that the magnetization changes 
isothermally at a temperature To. The actual 
TR (ITR) occurs if a dc-demagnetized mag- 

(2) 

net is heated (cooled) while H e3 ! t is kept con- 
stant. Changing the temperature has two ef- 
fects. On the one hand the saturation mag- 
netization decreases with increasing tempera- 
ture and on the other hand the switching fields 
H a are changed. Whereas the temperature de- 
pendence of the saturation magnetization can 
be measured easily, it is impossible to mea- 
sure it for the switching fields. Fortunately, 
one can solve eq. CU) for H s with H^J t = H c 
and (M) 1 = [23]. H C (T) and M S (T) are 
taken from measurements . If the influence of 
the weak grains is negligible, H S (T) is related 
to Hc(T) and Ms(T) by the simple formula 
HoH s (T) = fi H c (T)+nM s (T) 23 . Further- 
more we adopt the temperature dependence of 
H a for arbitrary values of the switching fields, 
i.e 



H S [T) = H S (T ) x (H S (T)/H S (T )) . 



(28) 



The calculation of the magnetization (M} 2 at 
an enhanced temperature T is similar to the 
calculation of the magnetization (M) 2 by help 
of eq. ()28j) . whereby the temperature induced 
deformation of the switching field distribution 
has to be regarded. Due to the assumption on 
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H S (T) and the normalization of the switching 
field distribution follows that a s has the same 
temperature dependence like H S (T). The mag- 
netization of a weak grain with given H s at 
the temperature T results from eq. (|14|) with 
the appropriate H S {T), Ms(T), and (M)^ in- 
serted. For the hard-grain contribution we have 
to calculate the probabilities w^, u>j|, w^, and 
!i)n, however for the increased (lowered in case 
of ITR) temperature T. We find for the mag- 
netization of the hard-grain fraction with given 

Mj(H s ) = M^{w n T (H s )-w n T (H s ) 

+w u T (H s )- Wll T (H s )y (29) 

It is obvious that the probability pi(H s ) 
remains unchanged, since the grains were 
switched downwards at the initial temperature 
To. Therefore we have 

w n T (H s ) =(l- Pl (H s ^(l-p^(H s ^ (30) 

™n T ( H s) = (l-pi(H s ))pT(H s ), (31) 
w n T (H s ) = Pl {H s )ql{H s ), (32) 
Wll T (H s ) = Pl (H s )(l-qJ(H s )). (33) 



The probabilities g<f (if s ) and p 2 {H s ) become 
temperature dependent both directly due to 
H S (T), M S (T), and <r/(T) and indirectly via 
(M) 2 . The temperature dependence of (Jf(T) 
was assumed to be that of Ms, since the fluc- 
tuations are caused by the inhomogeneities of 
the magnetization. Finally the TR or the ITR 
can be calculated from the following implicit 
self-consistent equation 

nM s 



{M) T 2 



dH s g(H s )(M) w /(H s ) 



CO 

+ J dH s g(H s )(M) 2 T (H s ). (34) 

nM s 

3.6 Influence of approximations 

The theory teaches us, which processes are 
necessary to describe different parts of the TR 



curve qualitatively right. In order to describe 
a residual remanence beyond Th c , Hs(T) has 
to be determined self-consistently from Hq{T) 
and M s (T). In Fig. l2"2l we study different 
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Figure 22: The TR for different levels of ap- 
proximation. The solid line gives the results 
of the presented theory. The dot-dashed line 
results from neglecting the weak grain frac- 
tion and for the calculation of the dashed line, 
both the weak grains and p 2 , i.e. the proba- 



bility, that a grain which resisted H^]} t 



will be 



(2) 

switched while heating at H^J t 
glected. 



0, were ne- 



approximation levels. The influence of weak 
grains is small as long as the temperature is 
below T max . For T > T max the neglect of 
the weak grains would result in a qualitative 
different behaviour, since we get a sign change 
of the TR before it finally vanishes, as shown 
in Fig. |2U in dashed-dotted manner. The 
reason for that peculiarity is made clear in Fig. 
12*31 where we have plotted the temperature 
dependence of the probabilities w^, w^, w^, 
and wn regarding the hard grain fraction. It 
is obvious that at first increases and when 
it drops down wu rises up, thus explaining the 
sign change. This behaviour is not observed 
if the probability p 2 , i.e. that a grain which 
resisted H^ x \ will be switched by H^ t , is 
neglected, as it was done in former theories, 
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but then an overestimation of the TR occurs, 
as visible in Fig. 1221 For the description of 
the TR at higher temperatures it is necessary 
to take into account, that above T max even 
the switching fields of the hardest grains are 
lowered enough to be switched by the still 
existing stray fields. The latter do not vanish 
until the magnetization breaks down. It is 
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Figure 23: The probabilities w^ T , w^ J ', w^ T , 

in dependence on 
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barium ferrite. Taking into account this pro- 
cess is in principle not diffult, since it simply 
scales the SFD of the switched particles. This 
would introduce at least one additional param- 
eter. We abstained from doing so by restrict- 
ing the theory to the first TR run only. Taking 
into account all probabilities and weak grains 
the model describes the TR dependence on the 
demagnetization factor N of the sample and 
the initial temperature To well, as was shown 
in Ref. |2H] for the TR and in Ref. [22] for the 
ITR. In the following we concentrate on the in- 
vestigation of the influence of internal model 
parameters. 

3.7 Comparison to the experiment 

The limitations of the presented model, i.e. 
the assumed perfect texture and single phase 
consistency, are best fulfilled in well aligned 
SmCo^- and barium ferrite magnets. There 
are three unknown parameters of the theory 
n, Of and o s which have to be determined by 
adjusting the theory to the experiment. Fig. 

0.4 



T 



T 



T 



T 



T 



and wii (cf. eqs. 

the temperature. The dashed line shows the 
ratio of the volume of weak grains to the total 
volume. 

worth to point out, that weak grains are not 
necessary to get effects of more than 50 %. 
Such a large TR one gets due to the feedback 
of the magnetization via the internal demag- 
netizing fields. In some sense this models 
the well known avalanche effects observed 
during hysteresis measurement, which should 
also occur to a certain extent if the sample 
remagnetizes. 

In the early works on the TR a further point 
was stressed jHlEHEl: The field dependence 
of switching field. This is a crucial point, if one 
attempts to describe repeating experiments or 
the virgin curve, but can be neglected in the 
first TR experiment, if H^) t is large enough 
to make most of the down-switched particles 

hard, as it seems to be the case in SmCo§ and [21 shows the best fit results for the SmCo§ 




700 800 
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Figure 24: Comparison of the best-fits for the 
SmCoz sample from Fig. [HJwith Hg(T) calcu- 
lated self-consistently (solid line) and approxi- 
mately (broken line) from Hc{T) and M${T). 
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sample with N=0.1 used in Fig. |H] where the 
(smoothed) measured temperature dependence 
of Hq and M5 (shown in Fig. 0) was used 
as input. The best agreement is achieved, 
if Hg(T) is calculated self-consistently from 
eq. (|18|). Although up to the maximum 
the difference is negligible the approximated 
formula fails at higher temperatures where 
most of the grains are weak. Furthermore the 
influence of the form of the m-dependence of 
<7/ was checked. It comes out that the initial 
slope is a little bit increased and the maximum 
becomes a little bit lower and broader if the 
function <r/(m) is chosen wider. Unfortunately 
these changes are to less in comparison with 
the inaccuracies of the present experimental 
data that a definite decision can not be made. 
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Figure 25: Calculated TR between To and two 
reference temperatures T re j in dependence on 
a small steady field applied while heating (solid 
lines) . The circles and squares are experimental 
points taken from Fig. |HJ 

Figs. |H1 and [23 relate the theory to the ex- 
periment. For that the parameters were deter- 
mined by fitting the standard experiment, i.e. 

(2) 

with H^ xt = 0. The graphs for the other values 

(2) 

of H ext 7^ were then calculated (not fitted!) 
with the same parameter set. The agreement 
is fairly well. Of course, it can be made per- 
fect by fitting the parameter for every curve 



separately. The variation of the parameters, 
given in Tab. El shows that the minima in the 
least-square sum are rather flat and that the fit 
data have to be taken as rough estimates. From 



H ext /mT 


AioO"s/T 


/i cr//T 


n 


VACOMAX 170 N=0.1 


+40 


2.526 


0.360 


0.299 





1.503 


0.467 


0.349 


-40 


2.324 


0.324 


0.256 


-80 


1.709 


0.417 


0.298 


-120 


2.200 


0.407 


0.268 


-240 


1.639 


0.421 


0.427 


VACOMAX 200 N=0.53 





1.60 


0.91 


0.56 



Table 2: The best-fit values for a VACOMAX 
170 sample, determined from the experimental 
curves given in Fig. |H]and for VACOMAX 200 
sample [16] . 

comparing the curves with other samples of the 
same material and by calculating the N- and 
To-dependence with the values given above we 
found the fit parameters belonging to H ex t = 
to be the most reliable ones. In Fig. [5^1 the 
theoretical curves are calculated in analogy to 
the measurements presented in Fig. 1171 From 
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Figure 26: Calculation of the influence of a 
small external field variation. Solid line: Best- 
fit TR curve for the same sample as in Fig. |SJ 
Dotted Line: Calculation for fj>oH ex t = ±5mT. 
Broken Line: Calculation for a for ^H ex t = 
zblOmT. 
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the comparison of both figures the qualitative 
agreement is ocularly, but the calculated sus- 
ceptibility is higher than the measured, what 
may be caused by the model assumption of 
ideal mobility of the domain walls within the 
weak grains. A further reason is the existence 
of a small volume fraction of lower coercivity 
components, which results in a small dip of the 
demagnetization curve, i.e. in a second broad 
flat peak of the H s distribution, which is not al- 
tered with temperature |2l]. That weak grains 
are responsible for the susceptibility enhance- 
ment becomes clear if one calculates the num- 
ber of weak grains in dependence on the tem- 
perature (dotted line in Fig. fTH|) . For the ITR 
in barium ferrite multi-domain grains are of less 
importance, since upon cooling the coercivity 
is not reduced below nMs if n is small. Thus 
most of the grains remain in SDS. Again we 
find an excellent agreement of the fitted curve 
with the experimental values, as may be seen 
from Fig. EJ From Tab. El The agreement is 



sample 


N 


a s /T 


a s /T 


n 


Bl 


0.1 


.124 


.128 


.156 


B2 


0.1 


.124 


.129 


.156 


B3 


0.3 


.120 


.106 


.153 


B4 


0.7 


.131 


.086 


.067 


B5 


0.9 


.102 


.113 


.049 



Table 3: Best fit values for the barium ferrite 
samples from Fig. E1(B1) and Fig. [7|(B2-B5). 

much better for the elongated samples, whereas 
for the prolate discs internal deviations from 
the assumed uniform magnetization seem to be 
more important |32| . In Fig. EZIthe calculated 
field dependence of the maximum TR for sam- 
ple Bl is shown, where the fit values from Tab. 
121 have been used. 

The influence of the intrinsic model parame- 
ters a s , <jf and n on the TR is studied in Fig. 
1281 for model parameters taken from Tab. [2] 
The most interesting features are the strong 
dependence on the field fluctuation width aj 
and on the internal demagnetization factor n. 
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Figure 27: Comparison of the calculated (solid 
line) and measured (diamonds) dependence of 
the ITR in barium ferrite on a small steady 
field applied while cooling from 300 K down to 
5 K. 

For small internal demagnetization factors, i.e. 
for elongated grains, the TR is small, whereas 
effects up to 100 % are possible for platelet- 
shaped grains. This corresponds well with the 
simulation results of Ref . [221 • The shift of the 
TR peak to higher temperatures as well as the 
sharpening with decreasing fluctuation width 
can be understood easily from the fact that 
the switching field of a volume fraction with 
given H s has to be reduced more by temper- 
ature before it contributes to the TR. Taking 
into account the different temperature depen- 
dence of He the same may be said for the ITR 
in barium ferrite, as seen in Fig. 1291 

4 Discussion 

Summarizing the experimental facts we find 
that TR or ITR is observed in all technical 
relevant permanent magnetic materials. Re- 
garding the prerequisites for TR and ITR 
we found that a necessary condition for the 
effect is the preparation of an asymmetrical 
demagnetized initial state, where the average 
switching field of positively and negatively 
magnetized grains is different. A further 
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Figure 28: Calculation of the influence of ay 
and n on the TR. The other parameters are 
the best-fit values from Tab. [21 



condition is the existence of a sufficient broad 
distribution of field fluctuations. Further- 
more the reduction of the coercivity with 
temperature has to be large enough, that the 
interaction fields may overcome the switching 
fields. The sign of the temperature coefficient 
of coercivity determines whether normal TR 
or ITR happens. Some features are common 
to all the four considered groups of perma- 
nent magnets (i) the strong dependence on 
the demagnetization factor ./V of the sample 
(ii) the extreme sensibility against a small 
superimposed steady-field H' (iii) the strong 
dependence on the initial temperature To (iv) 
the increase of the susceptibility after a TR 
cycle and (v) the reproducibility of TR in 
"repeating experiments". Otherwise, looking 
in detail, every permanent magnetic material 
shows some individual characteristics. For 
the SmCo^ magnets it is mainly the survival 
of a TR effect beyond the temperature Th c , 
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Figure 29: Calculation of the influence of n and 
df on the ITR. For the other parameters the 
best-fit values for sample Bl from Tab. El were 
used. 

due to the big difference Tc — Th c - For the 
NdFeB magnets it is the very small switching 
field distribution and the decoupling of the 
magnetic active grains which reduce the effect. 
The peculiarities of the barium ferrite mag- 
nets are caused by the different temperature 
dependence of the coercivity. The Sm^Coyj- 
type magnets differ from all the others by 
the underlying coercivity mechanism of the 
main magnetic phase. Regarding the latter 
materials, it is undoubtedly that the samples, 
where measurements were published, exhibit 
TR with all the main features. Otherwise the 
magnitude of the effect is so small, that it is at 
present impossible to decide, whether pinning 
controlled magnets really show TR. Moreover, 
also within a given material group, samples 
produced with different technologies may differ 
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considerably in the amount of TR. 

Nevertheless most of the observed features 
can be well understood in the framework of 
the above presented theory, at least qualita- 
tively, if it is possible to estimate how the spe- 
cial characteristics of a material are reflected in 
the model parameters. The latter are mainly 
the internal demagnetization factor n and the 
width of the field fluctuation cry. The param- 
eter a s results from the assumption of a gaus- 
sian SFD, which was introduced for simplicity. 
This simplification is of less importance if one 
calculates the TR, since it is very sensitive to 
field fluctuations but less sensitive to the SFD. 
For the hysteresis curve the situation is vice 
versa. Thus a single gaussian may not be accu- 
rate enough to fit the hysteresis curve in good 
quality. In principle this situation may be im- 
proved by calculating a more realistic SFD from 
the measured hysteresis curve. Regarding the 
calculation of the TR this would result in sub- 
stituting the erf-function by the adequate prim- 
itive of the SFD. We do not think that this is 
the main point, instead we regard the simple 
version of the "inclusion approximation" ( for 
details see [23]) as the most problematic ap- 
proximation, since it describes the retroaction 
of the different environments onto the different 
grains by only one parameter, the internal de- 
magnetization factor n. An argument in favour 
of this approximation is the fact, that it yields 
exactly the same relation between the coerciv- 
ity and the related (averaged) switching field 
as introduced both empirically [31] and theo- 
retically [SH]. Through this paper we used the 
term "grain" for the magnetic bistable units 
for simplicity. Indeed, these magnetic units 
may be different from the crystallites. Espe- 
cially if there are strong interactions between 
the grains, it may be that some grains agglom- 
erate and switch coherently. For instance the 
small internal demagnetization factor n in the 
above cited barium ferrite samples shows that 
the magnetic units are elongated, so it seems 
possible, that short chains of crystallites play 
the role of the magnetic units. The other ap- 



proximations are of less importance. Of course, 
the magnetic field fluctuates in three dimen- 
sions regardless of the assumed ideal texture, 
but, in the case of strong anisotropy it is merely 
the z-component of the field which switches the 
grains, and these fluctuations will be Gaus- 
sian again and are therefore contained in the 
model. The latter is supported by the TR and 
ITR experiments at isotropic sintered samples 
of SmCos and barium ferrite. Both the re- 
manence and the z-component of the fluctu- 
ating fields are reduced by a factor two with 
respect to the aligned samples, This yields ap- 
proximately the observed factor four for the re- 
duction of the TR and ITR resp.. The theory 
works best for the well textured SmCo^- and 
barium ferrite samples. For NdFeB, albeit 
curves may be fitted as well, the meaning of 
the fit parameters is not as clear. The existence 
of nonmagnetic phases will reduce the effective 
magnetization, i.e. the magnetic units may be 
thought ctS Si composition of a magnetic nucleus 
surrounded by a small nonmagnetic layer |19j . 
Thus the mean field can vary between ±cMs, 
with c being the content of magnetic phases, 
whereas in the switching conditions, it must be 
taken into account, that the switching unit is 
the magnetic nucleus only. Furthermore, the 
normalization of the SFD has to be changed 
to c. Without any calculation it may be un- 
derstood, that the TR in NdFeB is reduced 
in comparison with SmCo^. One reason is the 
smaller switching field distribution, which re- 
duces the asymmetry in the dc-demagnetized 
initial state. Furthermore, the nonmagnetic 
layer increases the distance between the mag- 
netic grains and therefore reduces the diver- 
gence of the magnetization, what results in 
smaller field fluctuations. If one does not as- 
cribe the observed TR in Sm,2Coyj to imper- 
fections, one has to take into account the cell 
structure yielding the volume pinning mecha- 
nism. If a grain of such a magnet is switched, 
the domain wall moves by successively switch- 
ing these cells, i.e. as switching units one has 
to regard now the cells, not the grains. Thus, 
also "single crystals" of Sni2Con should ex- 
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hibit TR. This was indeed observed [2U]. If 
the precipitation structure is perfect, all the 
cells will exhibit the same switching field, re- 
sulting in a very small switching field distribu- 
tion, what reduces the TR considerably. The 
switching of adjacent cells will be very corre- 
lated, resulting in a movement of the domain 
wall. Furthermore the exchange coupling of 
neighboring cells will surely dominate over the 
dipole-dipole interaction due to the smallness 
of the cells. Although this may be partly taken 
into account by allowing for internal demag- 
netization factors n bigger than one, the aver- 
aging over the neighboring grains as was done 
by the inclusion approximation becomes very 
questionable. Thus, the presented theory may 
be not very suited for the case of Sni2Coij. 

Regarding the data fitting we learned 
from Fig. [22] that both weak grains and 
down-switching of hard grains for higher tem- 
peratures has to be taken into account. This 
increased the computational effort consider- 
ably, due to the self-consistency requirements 
for H S (T) and (M). The model parameters 
determined from the fitting procedure have to 
be considered as estimates. To increase their 
reliability improvements of both the theory, 
as outlined above, and the experimental data 
are necessary. Regarding the experiment, both 
labor experience and the presented theory 
show that the most critical point in a TR ex- 
periment is the sample demagnetization factor. 
Of course, a closed circuit measurement would 
be highly desirable, as was shown in the very 
first reported TR experiment. This is hard 
to realize, especially for high-coercive samples 
and at low temperatures. Therefore rotational 
ellipsoidal samples with the symmetry axis 
along the texture axis should be used to secure 
the homogeneity of the mean internal field. A 
second point is a very accurate measurement of 
Hc(T) and Ms(T). Especially the coercivity 
should be measured before and after a TR 
experiment, where measurements beyond T max 
should be the last measurements at all, since 
they may result in irreversible changes within 
the material. Furthermore the field depen- 



dence of the irreversible susceptibility should 
be determined for several temperatures, since 
this would allow to determine an estimate of 
the SFD and to proof the assumption that the 
temperature dependence of the switching fields 
of different grains are equal. Since the theory 
was extended to arbitrary starting points, 
small deviations from the dc-demagnetized 
state are no problem. On the contrary, dif- 
ferent small steady fields should be applied 
during the TR experiment to determine the 
suppression field H', as an experimental 
estimate of the field fluctuation width. 

For the theory of the TR and ITR it was in- 
evitable to develop a theory for the hysteresis 
and the recoil curves also. Whereas there ex- 
ist a multitude of micro-magnetic calculations 
trying to explain the coercivity by assuming 
special micro-structures which give rise to ei- 
ther pinning of domain walls or nucleation, our 
theory starts a little above this level, since we 
do not judge about the reason for a switch- 
ing field, but simply accept its existence. In 
this point the theory resembles the Preisach 
HU model (PM). A theory of the TR based 
on the classical model is not able to explain 
TR effects of more than 50%. Furthermore this 
model was shown to have the congruency and 
wiping-out property [37] , which is a severe dis- 
advantage, since in most of the interesting cases 
the measured magnetization behavior violates 
these conditions. The main disadvantage of the 
classical Preisach model is due to the fact that 
the Preisach distribution function is not depen- 
dent on the magnetization state [HfiUUHj . To de- 
scribe TR effects up to 100% it was necessary to 
take into account many-particle effects, i.e. the 
feedback of the neighborhood onto the grain, 
which result in a dependence of the form of the 
fluctuation field distribution on the magnetiza- 
tion state ^01- As shown above the TR model 
was further completed step by step by taking 
into account the contribution of multi-domain 
particles and by the inclusion approximation. 
Parallel in time and mainly in the context of 
recording media the classical PM was improved 
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also by taking into account a mean field shift 
|23 EDI in tne moving PM (MPM). Although 
at first glance the "moving term" H + k(M) 
looks very similar, it describes a rigid mean 
field shift which is the same for all magnetic 
units. Contrary, the "inclusion correction" 
H — n(oMs — (M)))describes the deviation of 
the grain magnetization from the magnetiza- 
tion of its neighborhood. To take reversible 
processes into account the MPM was modi- 
fied another time by allowing non-rectangular 
elementary hysteresis loops |41| of the mag- 
netic units, or by splitting the Preisach den- 
sity function, as it was done in jEJEHl (GPM). 
The common drawback of all the mentioned 
Preisach models is the neglect of the state and 
history dependence of the form of the interact- 
ing field distribution. This had been shown al- 
ready by one of the first computer simulations 
|44j . and also by experiments |45j . The variable 
variance Preisach model accounts partly for 
that observation If the improvements of 

the various mentioned Preisach models would 
be mounted together, it may be probably pos- 
sible to describe the TR in permanent magnets 
also on that basis, if reasonable assumptions 
regarding the temperature dependence of the 
model parameters are added. A first step in 
this direction was done in Ref. |47j . 

In conclusion it can be said that the inves- 
tigation of the TR and ITR may be used to 
characterize modern permanent magnet mate- 
rials. In connection with the presented model it 
provides an alternative tool for studying inter- 
action effects in relation to the temperature de- 
pendence of the saturation magnetization and 
switching field distribution. 
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